Antimony arsenide: Chemical ordering in the compound SbAs 
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The semimetallic Group V elements display a wealth of correlated electron phenomena due to a 
small indirect band overlap that leads to relatively small, but equal, numbers of holes and electrons 
at the Fermi energy with high mobility. Their electronic bonding characteristics produce a unique 
crystal structure, the rhombohedral A7 structure, which accommodates lone pairs on each site. Here 
we show that the A7 structure can display chemical ordering of Sb and As, which were previously 
thought to mix randomly. Our structural characterization of the compound SbAs is performed by 
single-crystal and high-resolution synchrotron x-ray diffraction, and neutron and x-ray pair distri- 
bution function analysis. All least-squares refinements indicate ordering of Sb and As, resulting 
in a GeTe-type structure without inversion symmetry. High-temperature diffraction studies reveal 
an ordering transition around 550 K. Transport and infrared reflectivity measurements, along with 
first-principles calculations, confirm that SbAs is a semimetal, albeit with a direct band separa- 
tion larger than that of Sb or As. Because even subtle substitutions in the semimetals, notably 
Bii-^SW, can open semiconducting energy gaps, a further investigation of the interplay between 
chemical ordering and electronic structure on the A7 lattice is warranted. 

PACS numbers: 71.20.Gj 81.05.Bx 72.15.Jf 



I. INTRODUCTION 

The Group V elements As, Sb, and Bi have a structure 
all their own, unshared by any other elements. This "A7" 
crystal structure in Figure [TJa) is a consequence of the 
unique electronic structure of these elements , wh ere s-p 
hybridization leads to formation of a lone pair.QEl As a re- 
sult, the A7 structure has rhombohedral crystal symme- 
try and lies in space group i?3m. It is only a small distor- 
tion removed from simple cubic symmetry, which can be 
experimentally accessed under applied pressureP^ The 
band structure that drives formation of the A7 phase also 
causes these elements to be the prototypical semimetals, 
with a small offset band overlap, small number of carriers 
compared to typical metals (10 -5 as many), high mobil- 
ity, and nearly equal concentrations of electron and hole 
carriers P 

Their unique electronic properties have made the 
semimetals a fascinating arena in condensed matter 
physics, permitting initial measurements of the quan- 
tum mechanical oscillatory Shubnikov-de Haas and de 
Haas van Alphen effects, Seebeck's discovery of the 
thermoelectric effect, and Hall's measurements of spin- 
dependent transport P Because the band overlap in 
semimetals is so delicate, they are tunable by dop- 
ing for thermoelectric applications^ and in topological 
insulators P 

Alloying the Group V elements themselves can pro- 
duce unexpected results. The Bi 1 _ x Sb 2 , solid solution is 
semimetallic on both ends but becomes semiconducting 
for 0.07 < x < 0.22 by opening a gap at the L point 
and remo ving overlap at the T point of the Brillouin 
zoneP^EDBi and As are not chemically miscible. 12 Previ- 
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FIG. 1. (Color online) Unit cell of the A7 structure (a) of 
pure As, Sb, and Bi (and black P at high pressure) with R3m 
symmetry. Ordering the two 3a sites on this lattice leads to 
the GeTe structure in (b), lowering symmetry to R3m, while 
retaining the same two 3a sites. 



ous work on Sbi-^As^ has found that thes e elem ents are 
miscible across the full composition rangeP^S] Ohyama 
studied the thermal conductivity of Sbi_ x As x and spec- 
ulated that an anomaly around x = 0.5 might arise from 
an ordered compound but lacked any structural charac- 
terization to substantiate this claim.E^ 

While no pure elements outside Group V can form the 
A7 structure, the isoelectronic IV- VI compound GeTe 
is closely related. It is simply an ordered arrangement 
of the same structure, shown in Figure [jjb), with space 
group R3m in contrast to R3m for A7. GeTe might make 
a promising ferroelectric were it not for large concen- 
trations of free carriers arising from Ge vacancies and a 
Fermi energy lying within the valence band! 16 " 18 ! Again, 
the GeTe structure is a consequence of s-p hybridization 
and lone pairs on both Ge and TeP^l No other pure com- 
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pounds are known to crystallize in the GeTe structure, 
as nearby IV- VI compositions form the GeS or rocksalt 
structures, which are only small displacements awayP3 

Distinguishing between GeTe-type chemical order- 
ing (i?3m) versus an A7 solid solution (i?3m) is very 
difficult — no new Bragg peaks arise upon lowering sym- 
metry, so careful analysis of high-quality scattering data 
is required. For that reason, and to make clear the case 
for chemical ordering, we investigated SbAs using single 
crystal diffraction, high-resolution synchrotron powder 
diffraction, and pair distribution function (PDF) refine- 
ments. All find preference for 80/20 ordering of Sb and 
As. The PDFs do not show any short As-As distances 
which would signal the onset of nanoscalc phase sepa- 
ration that is invisible to Bragg diffractionPH^H High- 
temperature Bragg diffraction finds the atomic mixing 
to become disordered above 300° C. 

Transport and reflectivity measurements confirm that 
SbAs is a semimetal with band overlap, in agreement 
with our first-principles calculations. While the Seebeck 
coefficient and plasma frequency of SbAs lie between 
those of Sb and As, the optical dielectric constant is 
found to be outside the end members, implying a widen- 
ing of the direct band gap below the Fermi energy, in- 
dicating complex band structure changes. In the case of 
BiSb, this band shifting leads to opening of a semicon- 
ducting energy gap. None is found so far in SbAs, but the 
effects of doping and annealing remain a topic of further 
investigation. 



II. METHODS 

SbAs single crystals were prepared from elemental Sb 
(99.99%) and As (99.999%). The powders were loaded 
into 9 mm-diameter quartz tubes and sealed under vac- 
uum. Tubes were heated at 10°C/min to 800°C, at 
which point the samples were molten. Samples were held 
at this temperature for 30 min and periodically flipped 
to homogenize, then water quenched to avoid an incon- 
gruent melting transition.^. Still sealed under vacuum, 
tubes were placed into a furnace preheated to 630° C to 
anneal for 60 h, then cooled to room temperature at 
10°C/min. Under these conditions SbAs crystallized into 
shiny, mirror-like crystals about 1 mm per side, typically 
with triangular facets and easily cleaved into plates along 
{001} planes. 

Single crystal X-ray diffraction data were collected up 
to 9 = 35.28° on a STOE 2T image plate diffractome- 
ter equipped with Mo-Ka radiation (A = 0.71073 A) 
at room temperature. Data reduction and integration 
absorption correction were performed using X-Area soft- 
ware provided by STOE. The crystal structure was solved 
using direct methods and refined by a least- squares re- 
finement using the SHELXTL suite of programs.^ All 
atomic displacement parameters were refined anisotrop- 
ically. A twin law (-1-1001000-1) was applied 
and refined to 50.2%. The final composition refined to 



Sbo.94As1.06- Crystallographic parameters are given in 
Table ID 

High-resolution synchrotron X-ray powder diffraction 
was performed at beamline 11-BM of the Advanced Pho- 
ton Source (APS), using 30 keV x-rays (A = 0.413284 
A) and crystals ground and sieved to 45 /im. High- 
temperature diffraction was performed at beamline 1-BM 
using 20 keV x-rays (A = 0.6128 A) and samples sealed 
under vacuum in quartz capillaries. Time-of-flight neu- 
tron powder diffraction data were collected at the NPDF 
instrument at Los Alamos National Laboratory. Rietveld 
refinement was performed using the EXPGUI frontencP 
for GSAS. 28 . Unit cells and Fourier maps are plotted 
using VESTAP 

High-momentum-transfer total scattering data were 
collected at APS beamline 11-ID-B (90 keV, A = 0.13702 
A), and the aforementioned NPDF instrument (neu- 
tron time-of-flight). Extraction of the PDF was per- 
formed using PDFGetXd^and Q max = 25 A" 1 for x-ray 
data, and PDFGetlNP with Q max = 35 A" 1 for neutron 
data. Least-squares fits to the PDF were conducted with 
PDFguiP 

Resistivity measurements were performed in 4-point 
geometry using a Quantum Design PPMS. The Seebeck 
coefficient (thermopower) of a polycrystalline SbAs in- 
got was measured under helium atmosphere using an 
ULVAC-RIKO ZEM-3 system. Samples used for infrared 
(IR) reflectivity were annealed at 630° C for 30 h and 
cooled to room temperature in 60 h, resulting in large 
ingots. A flat surface of the sample was alumina pol- 
ished and washed with ethanol. The reflectivity spec- 
trum was recorded as a function of wavenumber, in nearly 
normal incidence, in the spectral range 100-4000 cm -1 
with a Nicolet 6700 FTIR spectrometer equipped with a 
Spectra- Tech spectral reflectometer. 

Electronic band structures and densities of states 
(DOS) were calculated for the hexagonal unit cell of 
SbAs using density functional theory (DFT). We used 
a perfectly ordered arrangement of SbAs but allowed 
for relaxation. All c alcula tions used the projector- 
augmented wave methocP^iJ anc j the generalized gradi- 
ent approximation to exchange correlation, developed by 
Perdew-Burke-ErnzerhofpSl as implemented in the VASP 
codeJMMl a planewave energy cutoff of 400 eV was used 
and convergence was assumed when the energy differ- 
ence between subsequent self-consistent cycles was less 
than 10 -4 eV. Self-consistent calculations were done us- 
ing 12x12x6 Monkhorst-Pack k-point sampling^ and 
the DOS is obtained by using a finer k-mesh of 18 x 18x9. 
Scalar relativistic effects and spin-orbit interactions were 
included. Thermopower S calculations using the Boltz- 
mann transport equation and constant relaxation time 
were performed with the BoltzTrap package written by 
Madsen and SinghP^ 
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TABLE I. Structural parameters obtained from room- 
temperature single-crystal refinement (full-matrix least- 
squares on F 2 ) of SbAs. fli = E||F | - |F C ||/E|F |, wR 2 = 
{■£[w(\Fo\ 2 -\F c \ 2 f]/Z[ w (\F f)]}^ 



Formula 
Formula Weight 
Crystal system 
Space group 
a = b 
c 

V,Z 

p 

Absorption coefficient 
F(000) 

8 max 

Reflections collected, unique 
Unique reflections 

Rint 

Number of parameters 
Goodness-of-fit on F 2 
Final R indices [I > 28(1)] 



Sbo.94As1.06 
193.86 g/mol 
Trigonal 
R3m 

4.0655(7) A 

10.889(3) A 

155.87(5) A 3 , 3 

6.196 g/cm 3 

28.748 mm" 1 

249 

34.62° 

762, 197 

197 

0.0203 

11 

1.311 

0.0117 Ri, 0.0287 wR 2 



III. RESULTS AND DISCUSSION 
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FIG. 2. (Color online) Rietveld refinement of SbAs using 
synchrotron X-ray radiation at beamline 11-BM of the APS. 
The full refinement range is given in (a), with the high-Q 
region enlarged to show detail. Ordering between Sb/As pro- 
duces a change in intensity in the {101} peak at Q= 1.88 
A -1 , highlighted and arrowed. This area is enlarged in (b), 
where the experimental and refined pattern is compared to 
the model using fully-ordered (dotted) and random solid so- 
lution (dashed). The difference between the two models (O-D 
difference) is shown at the bottom of (b). The largest inten- 
sity variation in the pattern lies on the {101} peak. 



Single-crystal diffraction provides the most observed 
reflections per refined parameter and so is the preferred 
method to determine whether SbAs forms an ordered 
compound or solid solution. The cell parameters and 
refinement details are summarized in Table [TJ while the 
atomic parameters are given in Table |LT} It must be noted 
that lowering symmetry from the A7 solid solution space 
group RZm to GeTe-type R3m does not result in the 
appearance of any new reflections. Only the peak in- 
tensities hold information about chemical ordering. We 
found that the cell contains two sites with different scat- 
tering density. The refinement indicated that an approx- 
imate 80/20 occupation of each site (see Table [ll| signif- 
icantly improves the fit versus a pure solid solution: re- 
finements with even site mixing in the R3m space group 
gave final R indices of R x = 0.0270 and wR 2 = 0.0754 
with a goodness-of-fit on F 2 of 1.527. This implies that 
SbAs has GeTe-type ordering with ^20% antisite disor- 
der. More complete ordering might be attainable by an- 
nealing below the ordering temperature which we discuss 
subsequently. 

Single crystal refinement indicates that SbAs is mostly 
ordered, but the two structures produce such similar 
diffraction patterns that multiple probes should be used 
to confirm correct occupancy. To that end, we performed 
synchrotron powder diffraction at beamline 11-BM of 
the APS, which provides unparalleled resolution and sig- 
nal/noise ratio for powder samples.^ The 11-BM data 



is shown in Figure [2] and the Rietveld fit is excellent. 
The 11-BM refinement comes to the same occupancies 
and Uij atomic displacement parameters as the single- 
crystal refinement, within the margins of experimental 
error (Table [TlJ) . 

Powder diffraction patterns allow us to visualize how 
much intensity is attributable to chemical ordering: the 
difference in the diffraction patterns for ordered and 
solid-solution SbAs models ("order-disorder difference") 
is plotted at the bottom of each pane in Figure [2] The 
difference was within the noise of conventional Cu-Ka 
diffraction data. The region with the largest order- 
disorder difference is highlighted in Figure ^ a) and mag- 
nified in Figure [2jb) . Here, the exceptional signal/noise 
ratio of 11-BM data provides distinction between the or- 
dered (dotted) and disordered (dashed) models in the 
{101} peak at Q = 1.88 A" 1 . From this view the 
increased experimental intensity versus the disordered 
model is clear, providing strong evidence for ordering. 

Since the diffraction pattern is the Fourier transform of 
the scattering potential and the experimental structure 
factors I -Fobs I can be extracted directly from the 11-BM 
data, a Fourier map can be constructed that contains 
the three-dimensional distribution of scattering density 
in the cell. 

The {100} slice of this map is attached to the side of 
the unit cell in Figure [3] to show the relationship between 
the Sb and As sites and their location on the map. The 
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TABLE II. Atomic parameters obtained from single-crystal and synchrotron powder refinement of SbAs. Atomic displacement 
parameters Uij are given in units of A 2 . For powder and PDF refinements, the total occupancy of each site was constrained to 
be 1, but the antisite fraction was not constrained to be identical for both sites. 
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FIG. 3. (Color online) The refined structure (a) from 11-BM 
X-ray data is shown with a Fourier map of observed structure 
factors Fobs on the {100} face. The slice is viewed normal to 
the plane in (b), with stronger scattering density evident on 
Sb versus As sites. 



{100} Fourier map viewed normal to the plane in Figure 
|3jb). Heavier scattering is evident on the Sb site [Z = 51 
versus Z = 33 for As). In a solid solution these two sites 
would have equal scattering density. 

The presence of nanoscale heterogeneity in SbAs can- 
not be entirely excluded by single crystal and powder re- 
finements, so we turned to PDF analysis. Local-structure 
PDF studies of the supposed solid solutions (In,Ga)As, 
Zn(Se,Te), and (Li,Na)AsSe2 have shown that materials 
with single-phase Bragg diffraction patterns can exhibit 
nanoscale clustering of the end members, evidenced by 
split nearest-neighbor bond distances in the PDFpSES 
Could there be nanoscale clustering of Sb-rich and As- 
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FIG. 4. (Color online) Least-squares refinements to the SbAs 
PDF collected at room temperature using (a,b) x-rays and 
(c) neutrons. In both cases, fits are excellent. The inset in 
(a) shows a two-phase mixture of Sb and As (dashed) that 
would indicate clustering. This is not present in the data. In 
the x-ray PDF, the calculated difference between ordered and 
solid-solution SbAs is larger than the noise. In the neutron 
case, due to similar scattering lengths, this difference is well 
below the noise threshold. 
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FIG. 5. (Color online) Sequential Rietveld refinements to 
high-temperature synchrotron x-ray diffraction data (beam- 
line 1-BM, APS) show disordering of Sb/As around 550 K, 
evidenced by the site occupancy in (a), while maintaining low 
goodness of fit R wp . In (b), the temperature dependence of 
the lattice parameters indicates a first-order transition around 
550 K. 



rich regions, as has been proposed by Levin, et al?^l For 
extensive nanoclustering of Sb and As we would see short 
bonds corresponding to As-As (2.52 A), plus long Sb-Sb 
bonds (2.91 A). The low-r region of the experimental x- 
ray PDF in Figure |4|a) shows a fit to the single-phase 
model from Rietveld refinement, plus a dashed line cor- 
responding to the nanoclustered model. None of the dis- 
tinct As-As or Sb-Sb bonds are present, eliminating the 
possibility of extensive nanoscale phase separation. 

A unit cell model can be least-squares refined using 
the PDF, just as was performed using single crystal and 
powder diffraction data. The x-ray PDF gives the same 
refined occupancy and atomic parameters as the single- 
crystal and powder refinements, providing a third check 
of the Sb/As ordering. The fit is shown in Figure ^b) 
and refined values are given in Table [TT] Neutron scat- 
tering does not provide sufficient contrast between Sb 
and As (5.57 and 6.58 fm scattering cross-sections, re- 
spectively) to resolve site ordering in the PDF or Bragg 
peaks, but higher r-space resolution in the neutron PDF 
(a result of higher usable Q m ax) further confirms the ab- 
sence of nanoscale phase separation by the good fit at 
low r in Figure |4]jc) . 

High-temperature diffraction data collected upon heat- 
ing a sample of SbAs at the APS beamline 1-BM refine 
to the same 80/20 site ordering as other techniques, and 
results of sequential refinements upon heating are shown 
in Figure [5] The site occupancy begins to deviate from 
the room-temperature value around 500 K. A value of 
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FIG. 6. (Color online) Resistivity (a) of SbAs shows impurity 
scattering at low temperature and linear temperature depen- 
dence, typical of semimetals. Resistivity of an ingot at higher 
temperatures (b) displays hysteresis, likely due to the phase 
transition around 550 K. Concurrent measurement of the See- 
beck coefficient 5* (c) shows p-type conductivity with a broad 
maximum around 500 K. 



0.5 in Figure [5ja) corresponds to even mixing (disorder) 
on both sites, and is marked with a dashed line. Fig- 
ure [5jb) shows a rapid change in the lattice parameters 
in the interval 500 K < T < 600 K, suggesting a first- 
order transition around 550 K. Annealing samples below 
this transition temperature may result in more complete 
chemical ordering. 



B. Transport measurements 

Resistivity measurements of SbAs crystals in the ab 
plane display metallic behavior with p — 80 /if2-cm at 
room temperature, shown in Figure[6ja) . The large resid- 
ual resistivity at low T is characteristic of impurity scat- 
tering. There is no region where the temperature coeffi- 
cient of resistivity is negative, as would be expected for 
a semiconducting region with finite E g . The behavior 
of SbAs can be contrasted with the Bii-^Sbz system, 
where 7-22% Sb substitution leads to an opening of E g 
up to 0.014 eV) 10 * 11 ! Saunders qualitatively suggested that 
SbAs does not become a semiconductor because the band 
overlap in As is greater than that in Bi, so a larger per- 
turbation from Sb addition would be required to shift the 
.L-point band enough to create a gapP^Figure[6|b) shows 
a rise in the resistivity above 500 K, coincident with the 
order-disorder transition observed by high-temperature 
diffraction in Figure [5] 

The measured thermopower of a polycrystalline SbAs 
ingot is shown in Figure |6^c). The positive (p-type) be- 
havior and magnitude are similar to that of pure Sb, 
which reaches a maximum of 30-55 uV/K at 400 K, de- 
pending on crystalline orientation!^! The Seebeck coef- 
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FIG. 7. (Color online) Experimental reflectivity spectrum (a) 
of SbAs and best-fit calculated reflectivity (solid line) using 
Equations [2] and [3] The Kramers-Kronig obtained Im(e) and 
the energy loss function Im(— 1/e) are shown in (b). The 
Drude fit (Equation [3| to Im(-T/e) gives the plasma fre- 
quency ujp and complex dielectric function e(uj). 



ficient of arsenic is about a factor of 4 smallerPS The 
decrease in S at high temperatures indicates an increase 
in n-type character, which could result from the effects 
of different mobilities for electrons and holes or changes 
in band overlap due to lattice expansion. Given the intri- 
cacies of the band structure near Ep in semimetals, we 
turn to IR measurements and DFT calculations to un- 
derstand how chemical ordering in SbAs may affect the 
band structure. 



C. Infrared reflectivity 

The experimental reflectivity spectrum of SbAs is 
shown in Figure [7] This spectrum is dominated by a 
structureless plasmon, yielding high reflectivity values in 
the low frequency range with a broad minimum around 
^1300 cm -1 , followed by almost constant reflectivity val- 
ues in the high frequency regime. The spectrum was an- 
alyzed by the Kramers-Kronig method to obtain the real 
£i and the imaginary e 2 parts of the complex dielectric 
function e{u). The Ei spectrum Im(e) and the energy 
loss function 



Im(l/e) =£ 2 /4 
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(i) 



are shown in Figure[7j The Im(e) spectrum is increased in 
the low frequency range denoting the contribution of free 
carriers, while the peak in Im(— 1/e) indicates a plasma 
frequency of ^1100 cm -1 . The latter represents the fre- 



TABLE III. Drude parameters derived from room- 
temperature reflectivity data of SbAs 

ui p (cm -1 ) Soo 7p (cm" 1 ) 7o (cm -1 ) 

1323 41.6 2044 1187 



quency of a longitudinal collective mode, when the entire 
carrier gas system is displaced relative to the fixed ions. 

The reflectivity R(u>) is expressed through the complex 
dielectric function as 



(2) 



We find that we can fit the measured reflectivity spec- 
trum of Figure |7^a) with a single-carrier double-damped 
Drude formula for the complex dielectric function^ 



~ i( Cfp ~ 7o)^ 
uj(uj + 170) 



(3) 



where is the optical dielectric constant associated 
with the bound electrons and uj p is the plasma frequency: 



Ne 2 



(4) 



where N is the free carrier concentration and m* is 
the carrier effective mass. In the typical Drude expres- 
sion for e(w), the free carrier damping factor j p = 1/r 
is considered constant throughout the entire frequency 
range. In the case of Equation [3j the carrier relax- 
ation time t is taken to be frequency dependent, giving 
a frequency dependent damping factor. Here j p repre- 
sents the linewidth of the plasma response centred at 
lj = uj p and 70 represents the linewidth of the absorption 
at uj = 0. The ratio 7 p /w p describes the motion of charge 
carriers: vibrational when the ratio is small, diffusive or 
incoherent when the ratio is large. Notice that Equation 
|3]reduces to the simple Drude expression when j p = 70^ 
Fitting the data in Figure [7] gives the parameters in Table 
|III| The single-carrier model fits our data well and also 
describes the IR reflectivity spectra of Bii ^Sb^ alloys 
in the semimetallic composition rangeP^^ This, how- 
ever, does not exclude a two-carrier system, i.e. electrons 
and holes, because a two-carrier Drude expression with 
plasma frequency given by Equation [5] is equivalent to 
that of a single carrier if the two carriers have the same 
relaxation time. 



Ne 

ml 



ml 



(5) 



For SbAs the plasma frequency uj p — 1323 cm 1 lies 
between the two end members: 1000 cm" - 1 for Sbj^and 
2419 or 2017 cm -1 for two different orientations of AsPl 
Using the values of Table |III| and Equation [4] we calcu- 
lated N/m* = 8.14 x 10 20 cm -3 for SbAs, which also lies 
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FIG. 8. (Color online) Electronic structure of SbAs (rela- 
tivistic and spin orbit interaction are included): (a) Band 
structure with a pseudogap along the YL direction; (b) pro- 
jected density of states (DOS) showing a finite value at Ep. 
Detailed band structures around Ef are shown below for (c) 
SbAs, (d) pure Sb, and (e) As. 



between the respective values for Sb and As (5.9 x 10 20 
and 8.9xl0 21 cm~ 3 , respectively) PUSH 

Meanwhile, the optical dielectric constant for SbAs 
was found to be 41.6, which is les s tha n the values for Sb 
and As (80 and 50, respectively).^^] The polarizability 
of the valence electrons relative to both parent materials 
is decreased, which could be explained by the empirical 
Moss relation, oc E~ x , which implies that de- 
creases as the direct band gap E g increases. This relation 
holds in Bii-^Sba, alloys across the band-opening compo- 
sition region .EIS2I A similar direct gap opening opening at 
the L point of SbAs, albeit small and below Ep, is found 
in our DFT electronic structure calculations of SbAs. 

Our IR reflectivity measurements confirm that SbAs 
displays carrier concentrations typical of Sb and As, with 
band overlap typical of a semimetal. However, the di- 
rect band separation probed by seems to have opened 
wider than that of the end members. We conducted first- 
principles calculations to visualize how chemical ordering 
affects the band structure. 



and electron pockets if the chemical potential is moved, 
for example by doping. 

The projected DOS in Figure |8jb) clearly shows the 
pseudogap structure with finite DOS at the Fermi level 
(Ep). From the projected DOS, one can see that there 
is strong hybridization between Sb and As orbitals. Hy- 
bridization of Sb s and As s orbitals gives rise to two 
bands in the DOS from -13.5 eV to -6.5 eV, separated by 
a gap of ~1 eV. The similar picture also applies to the 
hybridization between Sb p and As p orbitals. However, 
the separation between the bands is small, giving rise to 
an overlap region, forming the pseudogap. 

Comparing the band structure of SbAs with those of 
pure Sb and As (Figures [8jc-e)) reveals common features 
over a large energy scale, but salient differences appear 
near the pseudogap region. The splitting of the valence 
and conduction bands at the L point in SbAs is much 
larger compared to the pure compounds, which results 
in a larger direct band separation probed by IR reflec- 
tivity The effect of SOI in SbAs is stronger than those 
in pure Sb and As due to the lowering of symmetry to 
non-centrosymmetric i?3m. Our structural refinements 
confirm the presence of anti-site disorder in the SbAs 
samples which we measured transport and reflectivity. 
Further investigation is required to determine how this 
disorder affects the pseudogap structure and whether it 
can be tuned by annealing or doping. 

Starting from the DFT electronic band structure we 
calculated the thermopower S. For the nominally un- 
doped system we find that the thermopower is posi- 
tive 22 //V/K, in good agreement with our experimen- 
tal value of 25/i/K. However the calculated S is nearly 
T-independent in the range 300 K < T < 550 K, ris- 
ing from 22 to 26 /iV/K, whereas the experimental value 
increases from 25 to 40 /i/K. The origin of this strong 
T dependence is not clear. The effect of the structural 
transition around 550 K on the electronic structure near 
the Fermi energy and consequent change in the chemical 
potential \i with temperature and hence S merits further 
theoretical study. 



IV. CONCLUSIONS 



D. Electronic structure calculations 

The calculated band structure and DOS of SbAs are 
shown in Figure p[a,b). Our calculations clearly show 
that SbAs is a semimetal with a pseudogap, in agreement 
with resistivity and IR experiments. The detailed band 
structure along KTL near the pseudogap is given in the 
Figure |8jc) . The valence band maximum (VBM) and 
conduction band minimum (CBM) occur at ^0.3rL and 
^0.9rL resulting in two electron and two hole pockets. 
The band structure has an interesting feature near the T 
point where the VBM has a local minimum and the CBM 
has a local maximum. These bands could produce hole 



We find that the compound SbAs forms as an ordered 
GeTe-type structure with 80/20 ordering upon cooling 
from just below the solidus temperature, as confirmed 
by single-crystal x-ray diffraction, high-resolution syn- 
chrotron diffraction, and PDF refinements. The low-r 
PDF data precludes the possibility of nanoscale phase 
separation of Sb and As. Transport measurements con- 
firm semimetallic behavior analogous to the end mem- 
bers, with the exception of a direct gap splitting that is 
suggested by IR reflectivity to be larger than that of Sb 
or As. First-principles calculations indeed find an open- 
ing of the direct gap around the L point due to subtle 
changes in the p orbital hybridization caused by lower 
symmetry. 



The chemical order-disorder transition is found to 
be around 550 K by Rietveld refinements to high- 
temperature synchrotron diffraction data. This raises the 
possibility that improved ordering may be attained by 
slow-cooling or annealing below this temperature. The 
maximum in thermal conductivity observed by Ohyama 
corresponds to ordered SbAs, so the ability to tune this 
behavior by control of chemical ordering arises. Further- 
more, a reinvestigation of the Sbi-^As^ phase space is 
warranted. The related systems Bii-^Sb^ and Asi- x P x 
may also deserve more detailed study. While a wealth of 
experimental data finds no miscibility gap in Bii-^Sb^, 
the behavior of Asi-^P^ is comparatively unknown — 



contradictory reports of the hypothetical compound AsP 
have been summarized by Karakaya and Thompson!^ 
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